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SUMMARY 


The optimal control model (OCM) of the human pilot is applied to the 
study of aircraft handling qualities. Attention is focused primarily on 
longitudinal tasks. The modeling technique differs from previous applications 
of the OCM in that considerable effort is expended in simplifying the pilot/ve- 
hicle analysis. After briefly reviewing the OCM, a technique for modeling the 
pilot controlling higher order systems is introduced. Following this , a simple 
criterion for determining the susceptability of an aircraft to pilot induced 
oscillations (PIO) is formulated. Finally, a model-based metric for pilot 
rating prediction is discussed. The resulting modeling procedure provides a 
relatively simple, yet unified approach to the study of a variety of handling 
qualities problems. 


INTRODUCTION 


The advent of modern, digital stability and control augmentation systems 
has created a renewed interest in the study of aircraft longitudinal handling 
qualities. This renewed interest is attributable to two factors: First, 

the higher order nature of the dynamics typically associated with digital 
control systems make analytical prediction of handling qualities difficult. 
Contemporary handling qualities specifications^ are written assuming "classical” 
aircraft characteristics, e.g. , in the longitudinal mode, the existence of 
distinct and dominant short-period dynamics is assumed. With modern systems , 
the short-period characteristics may be dramatically euLtered by feedback and 
the higher order control system dynamics, themselves, may dominate the vehicle 
handling qualities. Second, shortcomings in predictive techniques are made 
even more critical by the fact that severe handling qualities deficiencies often 
arise in practice which are directly attributable to the higher order nature 
of the digital control law implementation. An example of this is the ability 
of high frequency phase lags or time delays in the control system to sharply 
degrade aircraft handling qualities and to be a contributing factor to pilot 
induced oscillations. ^ 

In the research to be described, a pilot modeling technique for handling 
qualities research, discussed in Ref. 3, is utilized and extended to cover 
higher order systems. The characteristics of over thirty aircraft configur- 
ations are analyzed, primarily in the longitudinal mode. Particular emphasis 
is placed upon those configurations where control system dynamics and time 
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delays have been recognized as contributing factors to handling qualities 
dificiencies. The contribution of vehicle/cont )1 system dynamics to PIO 
tendencies are outlined and a metric for pilot rating prediction is discussed. 


BACKGROUND 


The pilot modeling technique as discussed in Ref. 3 forms the frsnnework 
for the research described here. This technique utilizes the optimal control 
model of the human pilot euid a novel method for the a priori selection of 
dominant OCM parameters (index of performance weighting coefficients and 
observation noise/signal ratios). A brief tutorial review of the procedure 
for selecting index of performance weighting coefficients is now presented. 
Consider the longitudinal tracking task of Fig. 1 in which the pilot is 
attempting to minimize pitch attitude deviations 0(t) in the presence of 
atmospheric disturbances. Ignore the dashed "internal attitude command" 
for the present. An acceptable index of performance for this task would be^ 

j " E|^^-jrV2(t)/efi + 6'2(t)/6§]dt| . 

where 6(t) is control rate. 


As discussed in Ref. 3, we assign an^arbitrary maximum allowable deviation 
to the time rate of change of the error, oCt), and denote it 0„. Now an 
effective time constant T can be introduced to define maximum Viewable devia- 
tions of the integral and derivatives of 0(t) as: 


and 


'M 


M 


'M 


specified but arbitrary 

; 




( 1 ) 


-be Justification xor using a single time constant to represent the ratio of 
the maximum value of a variable to that of its next highest derivative rests 
upon the system bandwidth iniplications which follow when Eqs. (l) are used 
in implementing the OCM. We will also assign a maximum allowable deviation 
to the time rate of change of the pilot's control, S(t), and denote it iji^. 
Similar to Sqs. (l) we write 
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^1 " » 

*> to be selected ; 

» 

\ “ *’ ( 2 ) 


The value of is not arbitrary, however, but is found using Eqs. (l) and (2) 
and the vehicle dynamics as follows: Let the pitch attitude dynamics of the 
the aircraft be given by 



= K- 


. 1 - 1 ^ n-2 

s +a + 

n-2 

n . n-1 , 
s +b , s + * 
n-1 


• • +a^s+a 
•bj^s+b 


(3) 


Then, as explained in Ref. 3, we write 

[1/t"‘^ + lb + . . . + Ibil + IbplTj . 

®M * 

+ . . . + lail + la^lu 

Thus, once T is known, 6j^ and (and, if needed, 0j{, etc.) can be determined 
immediately. Choosing T involves selecting a domain of 1/T: 1 /Ut<1/T'^U/t 
and then plotting J, the value of the OCM index of performance, vs 1/T. The 
operating point or "knee" of this curve determines T. The knee is defined as 
the point where 


9J _ ’^It=:t/U ~ 

91og{l/T) “ ^^5 logCU/x) - log(lA7) 

Here Hr is a constant, nominally unity, which can be used to reflect manipulator 
characteristics, much like an efficiency factor; t is the pilot’s time delay 
(nominally 0*2 secs). /l is the value of the index of performance which 

results when T=t/4. 

The ability of the OCM parameter selection technique to provide a pilot 
model which matches measured pilot describing fuentions, remnant power spec- 
tral densities and root mean square (RMG) performance measures was demonstrated 
in Ref. 3. In addition, the modeling technique >ras shown capabu.e of providing 
qualitative and quantitative handling qualities assessments. The method for 
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selecting observation noise/signal ratios for the OCM is discussed in Ref. 3 
and will not be dealt with here. 


Although E(i. C3) shows dynamics of arbitrary order, 
dynamics of Ref. 3 were of the form: 

e, , 

jis) = ^ ^ 

s(s + 2 c 0) s+(0 ) 
n n n 


all the pitch attitude 


( 6 ) 


When higher order dynamics are encountered, the method for selecting the 
operating point needs to be modified slightly. The large phase lags typically 
associated with the dynamics of vehicles with higher order dynamics need to be 
reflected in choosing the domain of 1/T to be used in Eq. (5). To accomplish 
this, a delay Tj^ is defined as the delay which accrues when the vehicle 
dynamics of Eq. (3) are represented as 



Kg(s+l/Tj^)e~'^DS 

2 2 
s(s + 2 ? U) s+0) ) 
n n n 


( 7 ) 


The parameters on the right hand side of Eq. (j) are found usin- a program to 
fit a linear transfer function model to the actual vehicle dynamics.^ Equation 
( 5 ) is modified by simply replacing t with t+Tjj. The resulting equation is 
interpreted graphically in Fig. 2. Calculating Tjj and including it in Eq. ( 5 ) 
constitutes the extension of the methods of Ref. 3 to higher order systems. It 
is important to emphasize that the actual higher order vehicle dynamics are used 
in the modeling procedure; Eq. (7) is employed only to select Xp which, in 
turn, determines the domain of 1/T used in finding the index of performance 
weighting coefficients. 


APPLICATION TO AIRCRAFT LONGITUDINAL HANDLING QUALITIES 


Pilot Induced Oscillations 

Table I lists the aircraft configurations which have been analyzed in this 
study. The designations in the column 1 ?led "Configuration" use notation 
fo\u..J in the corresponding references. *ne first sixteen deal with high 
performance fighter-type aircraft in tracking or landing approach conditions 
and are t ten from Refs. 2, 5 and 6. These configurations constitute the 
test cases for the majority of the assessments. The next four configurations 
are taken from Ref. 7 and represent pilot-in-the-loop simulations of a hover- 
ing helicopter. Configurations 21-25 are flight test results from Ref. 9 in 
which the Princeton University Variable Response Aircraft (VRA) was used to 
determine the effect of digital sampling rates and time delays on longitudinal 
handling qualities. The vehicle dynamics appropr'.ate for IO 5 kts airspeed were 
used in the modeling procedure. Tlie pilot ratings which were used were average 
values obtained from altitude tracking and approach and landing tasks (Fig. 3 
of Ref. 9 ). Finally, configurations 26-32 are taken from Ref. iO where a 
moving base simulator experiment on the NASA Ames Flight Simulator for Advanced 
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Aircraft (FSAA) was described which investigated a wings level-turn control 
mode lor air-to-ground weapons delivery. Note that unlike the previous twenty 
five configurations , these involve lateral-directional aircraft handling q,uali- 
ties. The effective vehicle dynamics for the lateral gunsight aiming task 
were parameterized by a damping ratio undamped natural frequency 

and a pure time delay l^:ie data for the so-called "fine" task were 

used. This task is explained in Ref. 10. 

As an example of the modeling results. Fig. 3 shows the longitudinal 
open-loop pilot/vehicle characteristics for three of the configurations 

used in Ref. 2. Here, the NASA Dryden F-8^digital fly-by-wire aircraft is 
considered with a rudimentary augmentation system ("Pitch Direct") and three 
transport time delays of 0.13 sec, 0.23 sec, and 0.33 sec, respectively. The 
predicted effect of the time delays is apparent in the reduced open-loop 
crossover frequencies lOg. This open-loop characteristic obviously has a 
deleterious effect on the closed loop Q/Qq transfer functions as shown in Fig. 

4 (9/ej,=Y Y£./(l+YpYj,)). This transfer function is important is assessing PIO 
susceptability. Although the task has been defined as pitch-attittude distur- 
bance regulation, attitude commands 0^ internally generated by the pilot would 
be employed in precise altitude regulation (dashed line in Fig. l). Note 

in Fig. 4, that as tp increases, |6/9el and decrease at all frequencies. 

Perfect command following, of course, implies 0/9^=l.O at all frequencies. 

In Fig. 4, j0/9^,|<l.O for all configurations when o)<3.0 rad/sec, and is 
particularly poor for the configuration with tj/=0.33 sec. It can be readily 
shown that open-loop crossover frequencies less than 3-4 rad/sec will invariably 
result in poor closed-loop attitude command- following characteristics. The 
simplest and most direct way for the pilot to attempt to improve this closed 
loop command- following performance is to increase by increasing his static 
gain. If the pilot attempts this for the F-8 configuration with xp=0.33 sec, 
a very lightly damped closed-loop oscillation occurs at (»)=3.3 rad/sec (see 
Fig. 4). This is identical to the PIO frequency shown in Ref. 2 for this 
configuration. 

Similar results are also obtained for configurations from Ref. 5. Figure 

5 compares a pair of open-loop transfer functions obtained using configurations 
"11" and "12" from Ref. 5 and applying the pilot modeling technique discussed 
above. Once again, the dramatic difference in the crossover frequencies 

is apparent. The effects of the pilot's attempting to improve the performance 
of configuration "12" by increasing his static gain by 10 dBare shown in Fig. 6. 
Once again, a lightly damped oscillatory mode is seen to appear. The simula- 
tions of Ref. 5 were intended to provided performance comparisons for confi- 
gurations which vrere flight tested and discussed in Ref. 8. The latter report 
included Pilot-Induced-Oscillation-Ratings (PIOR) obtained using the scale of 
Fig. 7. It is interesting to note that configuration "ll" received an average 
PIOR or 1 indicating a very satisfactory vehicle whereas configuration "12" 
received a marginal average rating of 2.7 indicating a vehicle with definite 
PIO tendencies. These experimental results are seen to corroborate the ana- 
lytical findings Just discussed. 

Next, consider two configurations from Ref. 6 denoted as "4-1" and "6-1". 
Figure 8 shows the YpY^. plots for these configurations. Configuration "4-1’’ 
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received a very satisfactory PIOR of 1 whereas configuration ”6-1" received 
a very poor PIOR of U. Indeed, configuration "6-1" produced a PIO in flight 
test with a frequency of approximately 3.75 rad/sed. Analytically increasing 
the pilot's static gain by U.75 dB (the limit for closed-loop stability) in 
the modeling-results for this configuration produced a closed-loop oscillation 
at approximately 3-5 rad/sec. This 4.75 dB increase would increase oj^ from 
around 1.5 rad/sec to only around 2.5 rad/sec as compared to a value of 4.5 
rad/sec for configuration "4-1". 

Figure 9 shows the predicted YpY(.'s for a pair of configurations from 
Ref. 9. The task was longitudinal control in approach and landing using the 
Princeton VRA. The variable of interest here was the amount of effective delay 
in the control system. In the first, an effective delay of 0.055 sec was 
employed, while in the second, 0.355 sec was used. Again, note the striking 
difference in crossover frequencies in the predicted pilot/vehicle dynamics. 

In the first case, Uc=3.4 rad/sec, while in the latter, a)g=0.55 rad/sec. 

Flight test of the first configuration showed no PIO tendencies, while those 
for the latter produced PI0's.° 

Finally, Fig. 10 shows the predicted Y Y 's for a pair of configurations 
from Ref. 10. In the first, the control system parameters were ^.,=1.4, t 0 jj= 

2.0 rad/sec and Tp=0.0 sec, while in the second, 5^=1. 4, u)^«15.0 rad/sec and 
Tp=0.49 sec. The difference is again evident. Simulation results indicated 
that the configuration with delay was definitely PIO prone and the one without 
delay was not. It is interesting to point out that the configuration without 
delay still received an average Cooper-Harper pilot rating of 6.5, even though 
it was not PIO prone. Thus, poor pilot ratings, per se, are not a necessary 
condition for PIO susceptability. 

In each of the cases above, we have made direct comparisons of vehicles 
which were found to be PIO prone with those which were not. This was done 
to emphasize the fact that the method proposed here is clearly discriminatory 
in predicting PIO susceptability. The simple criterion for exonerating a 
vehicle from PIO tendencies requires that the predicted pilot/vehicle cross- 
over frequencies associated with inner attitude-loops be greater than 3-4 
rad/sec. 


Cooper-Harper Ratings 

Figure 11 is a plot of the Cooper-Harper ratings which the thirty -one configura- 
tions from Table I received in simulation or flight test /s, the value of a 
proposed handling qualities metric defined as K^* ((t+Tp)/t)^*J. No ratings 
were reported in Ref. 10 for configuration 32 of Table I. Hence, only thirty 
one data points are shown in Fig. 11. The Kj can be interpreted as a "cali- 
bration parameter" which, when multiplied by ( (t+Td)/t)^* J, allows the reported 
pilot ratings from different tasks and data sources to coalesce as shown in 
Fig. 11. Note that we do not allow Kp to vary within the analysis of any 
particular task, regardless of configuration changes. Thus, the analysis of 
the six configurations from Ref. 2 used a single value of (call it K,). 

The analysis of the seven configurations from Refs. 2 and 8 used a single value 


k 
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of K- (call it Kp), etc. In all, six different values (each one corres- 
ponding to the six different symbols in Fig. ll) were used .0 generate Fig. 

11. With the exception of all the parameters of the metric are an intrin- 
sic part of the modeling procedure, and, as such, involve no guessvork on the 
part of the analyst. In order to determine K^, the analyst must have an 
actual pilot rating for one of the configurations tested for the task under 
study. If the analyst does not have such a rating available, Fig. 11 is still 
useful, since the curve is nearly linear from a pilot rating of about 2.0 to 
10.0, a range which covers 80 % fo the Cooper-Harper scale. Thus, relative 
rating changes may be able to be predicted using the linear portion of the 
curve. Note that , with the exception of one data point (Config. 19 from Ref. 
7 ), the scatter in the ratings in Fig. 11 is only about *1/2 a pilot rating. 


The inclusion of the factor (.Ct+Tj))/t)^ in the metric deserves a brief 
discussion. In previous research with the OCM, the value of J, alone, has 
been found to correlate well with pilot opinion rating. In many of the 
configurations studied here, however (those with Tp>0) , the value of J was 
not acceptable as a metric. In general, the "predicted" opinion rating 
increments were smaller than those reported in experiment. There appears 
to be a reason for this based upon pilot tracking performance. Namely, when 
the task is disturbance regulation involving relatively low-bandwidth turbu- 
lence, large time de]ays are not necessarily a harbinger of dramatic deteri- 
oration in tracking performance. This is analytically verified by considering 
the R^B tracking scores for ccnfigurstions 1 and 3 from Table I. Here, a 13 ^% 
increase in time delay between configurations 1 and 3 invloves a logco^, 
regression of nearly a decade. However, the predicted RMS pitch attitude score 
increases by only 36 % and the predicted RMS control-rate score actually de- 
creases. As we have attempted to point out here, however, the same cannot 
be said for discrete command following or abrupt maneuvers. In this case, Ug 
regression can have a significant impact on the ability of the closed-loop 
pilot/vehicle system to follow abrupt, internally generated commands. It 
certainly is not unreasonable to postulate that such short-term response 
characteristics (in addition to R^B characteristics) are reflected in pilot 
opinion rating. Indeed, recorded pilot comments support this notion (e.g.. 
Refs. 2 and 6). The inclusion of ( (t+tjj)/t)** in the metric appears to account 
for the influence of these delays on pilot opinion in a straightforward manner, 
employing an easil" identifiable parameter (rp). 


CONCLUSIONS 


The research summarized in this paper provides a unified approach to 
pilox/vehicle analysis, and in particular for 

1 ) modeling the pilot controlling higher order systems 

2 ) prediotiag the susceptability of aircraft to longitudinal PIO's 

3 ) predicting pilot ratings for tasks when one configuration rating 

is known, or predicting relative rating changes between configurations. 
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Although the majority of tasks studied dealt with longitudinal control, five 
lateral-directional configurations were successfully analyzed with no changes 
in the modeling technique. 
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Table I. Aircraft Configurations Analyzed 


No. 


Configuration 


Reference 


1 

F-8 "Pitch Direct" 0.13 

sec delay 


r 

C 

> 

2 

.23 





3 

.33 





1» 

"ISAS" 0.13 

sec delay 




5 

.23 





6 

.33 



\ 

r 

7 

"2D" 



5.8 

8 

"5 A" 





9 

"OA" 





10 

ffgff 





11 

"10" 





12 

"11" 




1 

13 

"12" 




i 

Ik 

"3-1" 



1 

6 

15 

••1,-1" 



I 

1 

16 

"6-1" 



i 


17 

"PH-28" 




7 

18 

"PH-29" 





19 

"FH-32" 



1 


20 

"PH-35" 



\ 


21 

Princeton VRA 0.055 sec 

delay 


( 

9 

22 

.135 





23 

.255 





2*4 

• 355 





25 

.1i55 



\ 

r 

2b 

FSAA V/ings-Level Turn C ui 

(lateral-directional) r \ / 

1 rad/sec) 1 

1.1* 15.0 

sec) 

0 

10 

! 

27 

1.1* 

2.0 

0 



a:/ 

Cvi 

2.0 

8,0 

0 



29 

0.7 

6.0 

0 



30 

0.5 

1*.5 

0 



31 

0.3 

1*.5 

0 



32 

1.1* 

1*.5 

0.1*9 

\ 
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Figure 1. A Pitch Attitude Regulation Task 



Figure 2. Selecting an "Effective Time Constant" T 




a 

CO rad/Mc 

Figure 3. Pilot/Vehicle Dynamics for Three Configurations 
From Bef. 2. 

~ FOR tq - 0.33 see AND 5 dB INCREMENT IN 



rad/sec 

’^Igure I*. ''losed-Locj .'haracteristies for '^ree Configurations 
From Ref. 2 . 
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PIO TENDENCY RATING SCALE 


DESCRIPTION 

NUMERICAL 

RATING 

NO TENDENCY FOR PILOT TO INDUCE 

1 

UNDESIRABLE MOTIONS. 


UNDESIRABLE MOTIONS TEND TO OCCUR 
WHEN PILOT INITIATES ABRUPT 
MANEUVERS OR ATTEMPTS TIGHT 
CONTROL. THESE MOTIONS CAN BE 
PREVENTED OR ELIMINATED BY PILOT 
TECHNIQUE. 

2 

UNDESIRABLE MOTIONS EASILY 
INDUCED WHEN PILOT INITIATES ABRUPT 
MANEUVERS OR ATTEMPTS TIGHT 
CONTROL. THESE MOTIONS CAN BE 
PREVENTED OR ELIMINATED BUT ONLY 
AT SACRIFICE TO TASK PERFORMANCE 
OR THROUGH CONSIDERABLE PILOT 
ATTENTION AND EFFORT. 

3 

OSCILLATIONS TEND TO DEVELOP 
WHEN PILOT INITIATES ABRUPT 
MANEUVERS OR ATTEMPTS TIGHT 
CONTROL. PILOT MUST REDUCE GA<N OR 
ABANDON TASK TO RECOVER. 

4 

DIVERGENT OSCILLATIONS TEND TO 
DEVELOP WHEN PILOT INITIATES 
ABRUPT MANEUVERS OR ATTEMPTS 
TIGHT CONTROL. PILOT MUST OPEN 
LOOP BY RELEASING OR FREEZING THE ^ 
STICK. 

5 

DISTURBANCE OR NORMAL PILOT 
CONTROL MAY CAUSE DIVERGENT 
OSCILLATION. PILOT MUST OPEN 
CONTROL LOOP BY RELEASING OR 
FREEZING THE STICK. 

6 


Figure f. The Pilot-Induced-Oscillation Rating Scale. 
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Figure 10. Pilot/Vehicle Dynamics for Two Configurations 
From Ref. 10. 


CONFIGURATION 



Figure 11. Cooper Harper Pilot Ratings vs a Proposed 
Model-Based Metric. 
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